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Abstract 

We review and summarize recent theoretical and experimental work on electron spin 
dynamics in quantum dots and related nanostructures due to hyperfine interaction with 
surrounding nuclear spins. This topic is of particular interest with respect to several 
proposals for quantum information processing in solid state systems. Specifically, we in- 
vestigate the hyperfine interaction of an electron spin confined in a quantum dot in an 
s-type conduction band with the nuclear spins in the dot. This interaction is proportional 
to the square modulus of the electron wave function at the location of each nucleus leading 
to an inhomogeneous coupling, i.e. nuclei in different locations are coupled with different 
strength. In the case of an initially fully polarized nuclear spin system an exact analyt- 
ical solution for the spin dynamics can be found. For not completely polarized nuclei, 
approximation-free results can only be obtained numerically in sufficiently small systems. 
We compare these exact results with findings from several approximation strategies. 

1 Introduction 

In the recent years an extraordinary and increasing interest in spin-dependent phenom- 
ena in semiconductors has developed in the solid state physics community [1, 2]. These 
research activities are often labeled by the keyword "spintronics" which summarizes the 
entire multitude of efforts towards using the electron spin rather than, or in combination 
with, its charge for information processing, or, even more ambitious, quantum information 
processing. In fact, in the recent years a series of proposals for implementing quantum 
computation in solid state systems using electron and/or nuclear spins have been put 
forward [2, 3, 4, 5, 6, 7, 8]. In order to use the electron spin as an information carrier, 
long spin decoherence times are desirable if not indispensable. A serious possible limi- 
tation of spin coherence in semiconductors is the hyperfine interaction with surrounding 
nuclear spins. In fact, in semiconductors isotopes carrying a nonzero magnetic moment 
are ubiquitous. The commercial use of semiconductor technology is so far grossly domi- 
nated by silicon applications. In this material the magnetic isotope ^^Si having a spin 1/2 
and a magnetic moment of —0.5553 nuclear magnetons occurs with a natural abundance 
of 4.7%, apart from the two spinless stable isotopes ^*Si and '^"Si. From these numbers, 
hyperfine interaction might not appear to be a particularly relevant issue. However, the 
systems presently mostly studied in the field of spinelectronics and solid-state quantum in- 
formation processing contain materials such as GaAs, (Ga,Al)As and InAs whose elements 
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consist entirely of spin-carrying isotopes with substantial magnetic moments. As it will be 
discussed in detail in the present review article, in such systems the hyperfine coupling of 
electron spins to nuclear spins can easily become an important interaction. In table 1 we 
summarize the natural abundances and magnetic properties of stable spin-carrying nuclei 
relevant to semiconductor systems. 

A principle way to avoid such hyperfine couplings to electron spins is to use isotopi- 
cally purified material containing only a strongly reduced amount of magnetic isotopes 
[9]. However, present technology allows isotope purification of typical semiconductor ma- 
terials only up to a few hundredths of a percent or more of unwanted isotopes remaining. 
This degree of purification might in general not be sufficient to meet the high precision 
demands on implementations of quantum information processing. Moreover, such isotopi- 
cally purified materials have often quite high prices such that they might not appear as a 
technologically viable option even if the demands on precision are lower. 

This review is organized as follows: In section 2 we summarize important experimental 
work on electron spin dynamics in semiconductor nanostructures which have motivated 
the mostly theoretical studies we are reviewing in this article. In the following section we 
describe the basic physics of the hyperfine interaction between electron and nuclear spins 
as it occurs in semiconductors. Section 4 is devoted to a detailed analysis of electron spin 
dynamics in quantum dots due to hyperfine interaction with nuclear spins. After speci- 
fying the details of the model in section 4.1, we discuss the ground state and elementary 
properties of the underlying Hamiltonian in section 4.1.1. There we also describe a full 
analytical solution for the electron spin dynamics which can be obtained in the case of 
an initially fully polarized nuclear spin system. In section 4.1.2 we outline how to solve 
for the eigenstates and energy levels of the system via the Bethe ansatz technique relying 
on the integrability of the hyperfine coupling Hamiltonian. In section 4.2 we introduce 
different types of initial states for the nuclear spin system. These different types of initial 
states lead significantly different electron spin dynamics which are described in section 
4.3. The results reported on there are based on exact diagonalizations of the Hamiltonian 
for sufficiently small systems. In section 4.4 we discuss the intimate connection between 
the decay of the electron spin and the generation of quantum entanglement between the 
electron spin and the nuclear spin system. Section 4.5 is devoted to the important ques- 
tion of spin dephasing in an emsemble of dots as opposed to decoherence of of a single 
electron spin. In section 4.6 we discuss several further approximate treatments of the elec- 
tron spin dynamics that have appeared in the recent literature and compare them with 
the (quasi-)approximation-free approaches described before. Further, mostly theoretical, 
work relevant to the issue of hyperfine interaction between an electron spin bound to a 
quantum dot and surrounding nuclear spins is summarized in section 4.7. We close with 
conclusions and an outlook in section 5. 

2 Electron spin dynamics in quantum dots 

The dynamics of electron spins confined in semiconductor nanostructures is a rich and 
very active field. Among the most notable developments are experiments by Kikkawa and 
Awschalom who demonstrated very remarkably large coherence times for electron spins in 
n-doped bulk GaAs [10, 11]. This time scale can exceed 100ns and represents the time, 
i.e. the transverse spin relaxation time of an ensemble of electrons. These experimental 
findings have generated a great deal of prospects in the fields of spin electronics and 
semiconductor quantum computing [1, 2]. 

A photoluminescence study of excitons localized in single GaAs quantum dots was done 
by Gammon et al. [12]. In a subsequent study of the magnetooptical spectra of individual 
localized excitons [13] the role of electron spin hyperfine interaction with nuclear spins 
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was investigated. Though these experiments dealt with single GaAs quantum dots (not 
ensembles of them), it was not the electron spin dynamics studied there, but spectra of 
individual localized excitons where the strong Coulomb interaction between an electron 
and a heavy hole is important. The relaxation lifetime of electron spins for an ensemble of 
CdSe quantum dots of very small diameter (20-80 A) was measured by Gupta et al. using 
a femtosecond-resolved Faraday rotation technique [14]. Again, due to the small size of 
the dot the Coulomb interaction between electron and hole is not neglegible. Therefore, it 
is highly probable that the observed short relaxation time (being of order of nanoseconds) 
is due to fast spin dynamics of the hole. The surface states could also contribute to the 
spin relaxation due to the small size of the dot. Finally Epstein et al. have analyzed the 
spin lifetime of photogenerated carriers in InAs quantum dots using the Hanle effect [15]. 

The above experiments [14, 15] were performed on ensembles of quantum dots, not 
single dots. Therefore the estimated spin relaxation times are Tg* time scales. The experi- 
ments which directly probe the single electron spin relaxation were done recently by Fuji- 
sawa et al. [16] and by Hanson et al. [17]. The non-equilibrium tunneling current through 
excited states in an AlGaAs/GaAs quantum dot was studied using a pulse-excitation tech- 
nique which measures the energy relaxation time from the excited state to the ground state. 
Very low spin-flip rates were observed which is consistent with the theoretical predictions 
[18, 19, 20]. The physical mechanisms for electron spin relaxation for delocalized states in- 
clude the interplay of spin-orbit coupling with impurity scattering and/or electron-phonon 
interaction, and the hyperfine interaction with surrounding nuclear spins. In recent the- 
orctic;aI studies Khaetskii and Nazarov have concluded that the first type of mechanisms 
is strongly suppressed for electrons localized in quantum dots [19, 20], see also Ref.[21]. 
As to the contribution to the electron spin decoherence due to the combined effect of the 
spin-orbit interaction and the spin-independent interaction with acoustic phonons, there 
is an indication that within the Markovian approximation, which is usually applicable to 
this problem, T2 time can be as long as Ti time. This is due to the fact that an additional 
contribution to the I/T2 rate originating from the fluctuations of the spin-orbit related 
magnetic field along the external magnetic field direction is proportional to the phonon 
density of states at zero frequency. Thus, this contribution is zero for acoustic phonons, 
see Ref.[22]. These results have motivated a whole variety of theoretical investigations on 
hyperfine interaction in quantum dots and related structures which will be reviewed in 
this article [23, 24, 25, 22, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40]. 

A scenario similar to an electron bound in a quantum dot is the case of a shallow 
phosphorus donor in a silicon crystal, Si:P. This case is essential for the solid state quantum 
computing proposal by Kane [5]. The donor electron is bound to the P atom in a large 
hydrogen-like orbit with a Bohr radius of about 30A. The nuclear spins interacting with 
the electron are the central '^^P and the surrounding ^^Si. The Ti time scale for energy 
relaxation of the electron spin was determined by Feher and Gere to be of order lO^s [41]. 
The transverse spin decay was investigated by Gordon and Bowers [42] using the spin-echo 
method deducing a time scale of order 500/US, see also [43]. A very recent spin-echo 
study by Tyryshkin et al. [44] on P donors in natural and isotopically purified Si has 
reported T2 times being significantly larger than the previous result. 

Finally we mention that hyperfine interaction between electron and nuclear spins is 
of course also investigated in higher- dimensional semiconductor nanostructures such as 
quantum wells. For recent work using optical NMR techniques we refer to [45, 46]. 

3 Hyperfine interaction in semiconductors 

Hyperfine interaction is the coupling of a nuclear magnetic moment to the magnetic field 
provided by the (orbital and spin) magnetic moment of electrons. The Hamiltonian de- 
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scribing this interaction in the as a lowest-order relativistic correction was clcrivcid in 1930 
by Fermi [47]. For an s-electron there is no orbital contribution, and the Hamiltonian 
reads [47, 48] 

n = ^i,Bi^i\Wi)?sT. (1) 

Here S is the spin of the electron, and V'(rr) is its wave function at the location f/ of 
the nucleus. This Hamiltonian couples the electron spin to the nuclear spin I with total 
spin quantum number / and magnetic moment which is represented by the operator 
M/ = Both spin operators are taken to be dimensionless, /xo = 47r • 10~^Vs/Am 

is the usual magnetic constant in SI units, and /Us is the Bohr magneton. The leading 
contribution to the hyperfine coupling for higher angular momenta of the electron looks 
different from (1) and is essentially given by the usual dipolar coupling between the nuclear 
magnetic moment and the (orbital and spin) magnetic moment of the electron. 

For the issue of spin coherence and interaction with nuclei in semiconductors, the most 
relevant case are electrons in s-type conduction bands. We therefore shall concentrate 
on this case where the hyperfine coupling is described by the Hamiltonian (1). Prom 
the viewpoint of nuclear magnetic resonance experiments this coupling is the origin of an 
increment in the position of resonance lines known as the Knight shift. From the viewpoint 
of electron spin resonance effects in solids this Hamiltonian describes the Overhauser field. 
As a general reference on both effects we refer the reader to the textbooks by Abragam 
[49] and by Slichter [50]. 

In a semiconductor crystal the electron wave function is a product of a Bloch amplitude 
u{r) and a modulating envelope function ^{r), i^if) = ^{f)u{r). We therefore can rewrite 
the Hamiltonian as [48] 

H=^^lB^^IV\^{rI)\'SI (2) 

with rj = \n{ri)\'^ . For free electrons the Bloch function is constant, |u(r)| = 1; in a realistic 
crystal |M(r)| has maxima at the lattice positions, i.e. the locations of the nuclei, leading 
to ?7 > 1. Values for rj in semiconductor systems can be estimated from electron spin 
resonance experiments [48]. For InSb Gueron [51] found rjin = 6.3 • 10'^, ?7sb — 10.9 • 10^; 
for GaAs Paget et al. [52] estimated ?7Ga = 2.7-10'^, ?7as = 4.5-10'^. From NMR experiments 
on ^^Si Shulman and Wyluda estimated a value of jjsi = 186 [53]. 



4 Electron spin dynamics in quantum dots due to hy- 
perfine interaction with nuclei 

4.1 Modeling hyperfine interaction in quantum dots 

We consider an electron confined in a semiconductor quantum dot in an s-type conduction 
band. We assume the electron to be in some orbital eigenstate according to the confining 
potential, e.g. the orbital ground state in the quantum dot. The remaining spin degree of 
freedom is coupled to an external magnetic field B with an electronic g- factor g, and to 
the spins of surrounding nuclei via the the hyperfine contact interaction described in the 
previous section. Thus the Hamiltonian reads 

H=9iiBSB + sY,^ih- (3) 

i 

Here the subscript i labels the nuclei, and the coupling constants Ai are given by (cf. (2)) 

Ai = Avo\^{n)\'' (4) 
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where no = l/i>o is the density of nuclei. Provided that the electronic envelope wave func- 
tion \E'(r) varies smoothly on the length scale given by ^^Jvq it is appropriate to replace 
the sum Ylii by an integral over space; then A = ^ . Ai up to small corrections to this 
approximation. The Hamiltonian (3) also describes the hyperfine interaction between nu- 
clear spins and the spin of an electron bound in a hydrogen-type orbit around a phosphorus 
donor in a silicon crystal. 

Let us now address the order of magnitude of the hyperfine interaction. A GaAs 
quantum dot with a volume of order lO^nm^ contains of order N = 10^ nuclei with a 
density no = 45.6nm~^. Taking into account the natural abundances of the three occurring 
isotopes ^^Ga, ^^Ga, ''^As one has an average nuclear magnetic moment of /x/ = 1.84/X]v 
With the values for 77 estimated in [52] this leads to an overall coupling constant A of 
order 10~^eV to 10~^eV. This is the strength of the hyperfine coupling acting on the 
electron spin in the presence of a fully polarized nuclear spin system. With the effective 
electron g-factor g = —0.44 for GaAs this corresponds to an effective magnetic field of 
order a few Tesla. For a completely mipolarized nuclear spin system the strength of the 
hyperfine field is fluc;tuating around zero with a typical value given by A/^/N where N 
is the number of nuclei in the dot effectively interacting with the electron spin. Note 
that for GaAs all nuclear magnetic moments are positive, leading to an antiferromagnetic 
sign for the hyperfine coupling. This is different from the situation in Si:P where the 
magnetic moment of "^^P is positive while the surrounding ^°Si have a negative nuclear 
magnetic moment, resulting in a dominantly ferromagnetic coupling to the electron spin. 
The natural abmidance of ^''Si leads to a density no = 2.3nm^'\ Assuming a Bohr radius 
of 30A for the hydrogen- like electron orbit one can estimate the number of nuclei effectively 
interacting with the electron spin to be of order a few hundred, and from (5) one finds 
values of \A\ of order 10~'^eV using the estimate for rj given in [53]. 

Another kind of interaction in the systems discussed above is the dipolar coupling 
between nuclear spins. This contribution to the Hamiltonian is obtained from its classical 
counterpart by expressing the magnetic moments in terms of nuclear spin operators, p,i = 
{^1/1)1. Thus the interaction between two nuclear spins labeled by 1 and 2 reads 



where r\2 is the distance vector between the nuclei. In GaAs ^^Ga has the largest nuclear 
magnetic moment with jii = 2.562/xjv. Assuming two of such nuclei being on nearest 
neighbour sites of the underlying zinc-blende lattice one finds ri2 = 0.24nm, and the energy 
scale of the dipolar interaction between these two nuclei has the value (Mo/47r)y^j/rf2 = 
7.6-10~^^eV. This is an upper bound for the dipolar coupling between neighboring nuclei in 
GaAs and sets the time scale (10~^ 10~^s) on which this interaction indirectly influences 
the quantum dynamics of the electron spin. As we shall see below, the largest time scales 
relevant for electron spin dynamics due hyperfine coupling axe of order 10~^s. Therefore, 
the time scale of the dipolar coupling is much larger than the time scale provided by 
the hyperfine interaction; analogous considerations can be made for the case of Si:P. In 
the following we shall therefore neglect the dipolar interaction unless stated otherwise. A 
contribution to the Hamiltonian which we also neglect is the coupling of the nuclear spins 
to the external magnetic field. This interaction is much smaller than the Zeeman coupling 
of the electron spin because of the smallness of the nuclear magneton compared with the 
Bohr magneton for electrons. 



W12 = — 



A^O /<-lM2 1 / 3(Jifi2)(j2n2) 

47r /1/2 r?2 V ''12 



-/1/2 



(6) 
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4.1.1 Ground state and elementary properties 

If the sign of the hyperfine interaction is ferromagnetic {A < 0, as realized by the coupHng 

to ^^Si nuclear spins) the ground state of the Hamiltonian (3) is (for B = 0) the fully 
spin-polarized multiplet with all (electron and nuclear) spins in parallel resulting in the 
maximum value of the total spin quantum number, J = N 1+1/2. For an antiferromagnetic 
sign the classical ground state (neglecting the operator nature of spins) has all nuclear spins 
in parallel and the electron spin pointing opposite to them. We therefore anticipate that 
the true quantum mechanical ground state will have (again for vanishing external magnetic 
field) the total spin quantum number J = NI — 1/2. This assumption is confirmed by our 
numerics. Fig. 1 shows the energy spectrum of the Hamiltonian 3 for zero magnetic field 
as a function of the total spin quantum number J in a system of iV = 13 nuclear spins of 
length I = 1/2. Clearly the ground state lies in the subspace of Jz = (13 — l)/2 = 6. A 
general state in this subspace can be written as 

\^)=aS-\it,V--1) + ^^init,V--1) (7) 

where i labels the nuclei, | fl", T • • • T) is the fully spin-polarized ground state with all spins 
pointing upward, and we have introduced the usual spin lowering operators for the electron 
spin of length 1/2 and the nuclear spins of length /. The stationary Schrodinger equation 
leads to the following system of equation for the amplitudes a and 



(8) 



iA,a+(^^{A-2A,) + '-^-E^/3, = (9) 

Here E is the energy eigenvalue, and we have reintroduced a finite Zeeman coupling e^. = 
qIj^bB. Prom the second of the above equations one finds 



E-{IiA-2A,)+ez)/2- 

The ground state energy for je^j <C A is of order E « —AI/2, thus the denominator of the 
r.h.s. of (10) is of order A. Since the couplings Ai are of order A/N we see that each /3i 
is smaller than a by a factor of order 1/N. Therefore |a| « 1 up to quantum corrections, 
and all coefficients (3i{E) are of order 1/A^. Thus, the corrections to the classical ground 
state {\a\ = 1) are of order 1/A^, |ap = 1 — J2i lAP ~ 1 — 1/N. In summary, for large 
systems {N ^ 1), the ground state is essentially given by a tensor product state with all 
nuclear spins in parallel and the electron spin pointing opposite to them. Excited states 
with the same total spin quantum number along the direction of the nuclear polarization 
are separated from the ground state by a substantial gap of the order of the coupling 
parameter A. 

These results were obtained recently by Khaetskii, Loss, and Glaznian [24] who studied 
the time evolution of the classical ground state under the quantum Hamiltonian (3). These 
investigations were carried out using the standard Laplace transform technique. A detailed 
account on the mathematical details of this approach has been given recently in [25] . The 
main findings arc the following. The classical ground state for [e^l <^ A remains constant 
in time up to quantum corrections of order 1/A^. Starting the time evolution at t = 
with \a{t = 0)1 = 1 the electron spin expressed in terms of its expectation value {S^{t)) 
undergoes coherent oscillations between {S^{t)) = —1/2 and (S'^(t)) = —1/2 -|- 0{1/N) 
with a period of T = A-kU/A over a time scale of order %N/A. This time scale is nothing but 



6 



the; characteristic period of precession of an individual nuclear spin in the field generated 
by the electron spin. At this time scale a different nuclear spin configuration is created, 
and because of the spatial variation of the hyperfine coupling constants inside the dot, this 
leads to a different random value of the nuclear field seen by the electron spin and thus 
to its decoherence. After a time interval of this order the oscillations fade out and the 
expectation value {S^{t)) decays to a value of order {S^{t — > oo)) = — 1/2 + 0{1/N), see 
Fig. 2. 

Thus the decaying part of the initial spin state has smallness 1/N which is due to 
a large gap ~ A seen by the electron spin through the hyperfine interaction for a fully 
polarized state. As a result, only a small portion ^ 1/A^ of the opposite (+1/2) state can 
be admixed. Moreover, the decay of the electron spin turns out to be nonexponential for 
any external magnetic field. In the case when the perturbative treatment is applicable (i.e. 
a large Zeeman field) it follows a power law with a leading term proportional to t~^/'^. 
The fact that the decay of the electron spin is not exponential can be easily understood. 
The exponential decay occurs when the correlation time of the randomly fluctuating field 
which causes the decoherence is short compared to the decoherence time, as a result the 
Markovian approximation can be applied. In our case the nonexponential behavior is a 
result of the fact that the correlation time for the nuclear magnetic field seen by the electron 
spin is itself determined by the flip)-flop processes since the internal nuclear dynamics is 
excluded. Thus the Markovian approximation is not valid, see for comparison section 4.6.2. 

A particular situation arises in the time evolution of the state (7) with initially \a{t = 
0)1 = 1 if an external Zeeman field is applied to the electron spin which approximately 
cancels the initial Overhauser field, i.e. « —AI in the above conventions [24]. Near this 
Zeeman field |ap averaged over time is 1/2, i.e. the up and down states of the electron spin 
are strongly coupled via the nuclei (see Fig. 3). In contrast, outside this resonance regime 
the value of jap is close to unity (with small 1/N corrections), i.e. {S^{t)) = 1/2 — \a\'^ is 
close to -1/2 at any time. The width of the resonance is ~ A/\/N, i.e. small compared to 
the initial gap AI. This abrupt change in the amplitude of oscillations of {S^'{t)) (when 
changing in a narrow interval around AI) can be used for an experimental detection of 
the fully polarized state. 

Although physically not particularly realistic, it is also instructive to study the model 
(3) in the case of all coupling constants being equal to each other, Ai = A/N for all i. 
Models of this type were studied recently by Eto [33] and by Semenov and Kim [34] , see 
also [25]. The technical advantage of this type of models lies in the fact that the square 
of the total nuclear spin Itot = J2i an additional conserved quantity, [H, I^gf] = 0. If 
Ai = A/N the Hamiltonian reads in the absence of an external magnetic field 

Since the total spin J = Itot + S can have values J = hot ±1/2, each value of the quantum 
number Itot corresponds to two energy levels given by 

£(/)-5^(±(/« + i)-i). (12) 

These energy levels are typically highly degenerate. For instance, if the nuclear spins are 
of length 1/2, the different values of Itot occur with a degeneracy of 

( ^ \( ^ 
\N/2 - Itot) \N/2 - Itot - 1 

As seen from Eq. (12), the spectrum is equidistant with a level spacing AE — A/2N. 
Therefore the time evolution of an arbitrary state is strictly periodic with a recurrence 
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time T = AnUN/A. However, in contrast to naive expectations even in this simple case 
there is some time dependence of {S^{t)) which cannot be described by a single frequency. 
Actually for a nuclear state with given 7^^^ the solution contains all the frequencies of the 
form AE = A{2Iioi + 1)/2N, where Itot are all the moduli of the total nuclear momentum 
which can have this projection 1^^^. 

Depending on the number and length of the nuclear spins, this can lead to shorter 
periodicities in the time evolution of {S{t)). For instance, for an odd number of half- 
integer nuclear spins 2Itot + 1 is even, and {S{t)) has a period of T = 2TrJiN/A 



4.1.2 Integrability 

If all nuclear spins are of length 1/2, the spin Hamiltonian (3) has the strong mathematical 
property of being intcgrablc and solvable via an appropriate version of the Bcthe ansatz. 
This fact was recognized first by Gaudin in a rather formal context [54]. Consider a system 
of N spins 1/2 and the following sequence of operators [54, 55] 

where Oi are Pauli matrices and the Zi are some arbitrary coupling parameters. Obviously, 
by fixing one spin i to be the central spin and adjusting the couplings to the other spins, 
the operator Ti^ assumes the form of the Hamiltonian (3). The operators commute 
with each other, 

[Wi,Wj]=0, (15) 

and fulfill the sum rule 

i 

Thus, N — \ oi these operators together with the square of the total spin form a set of N 

linearly independent commuting operators being bilinear in the individual spin operators. 
The coordinate Bethe ansatz diagonalizing simultaneously all Jii can be summarized as 
follows [54, 55] . Consider states of the form 

\w^,...,w^)^F{w^)---F(w„,)\ t ••• t) (17) 

where | t • • • t) is the spin-polarized state with all spins in parallel, and 

^(-) = E^ (18) 

i 

with (j~ = (jf — iaf . These states are eigenstates of all Hi if and only if the complex 
parameters wi . . . Wm fulfill the Bethe equations 

JV m 

V^— + V^— =0. (19) 

Wk- Zi j-" Wl-Wk 

For a given solution of these equations the corresponding eigenvalues read 

£i(«;i,...,«;J =^^— + ^^— . (20) 

J ^ W Zi ... Zi Zj 

The eigenstates (17) are states of highest weight in multiplets of the total spin with respect 
to its 2;-component. States lower in these multiplets can be obtained by applying the 
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lowering operator of the total spin, or by formally considering solutions to the Bethe 
equations (19) with some Wk being infinite. By counting the number of solutions of (19) 
it can be shown that the Bethe absatz produces all multiplets, i.e. all energy eigenvalues. 

The coordinate Bethe ansatz outlined above has been extended to the technique of 
the algebraic Bethe ansatz [56]. We also note that the equations (19) are a limiting 
case of the Bethe equations of the so-called Richardson model describing electron pairing 
in superconducting grains [57]. This issue has attracted considerable interest recently 
[58, 59, 60, 61], and an algebraic version of the Bethe ansatz has also been presented [60]. 

For practical purposes, the treatment of the Bethe equations (19) is still rather com- 
plicated, and explicit results, for instance for correlation functions, arc difficult to obtain 
[62]. In the remainder of this review we shall therefore concentrate on techniques other 
than the Bethe ansatz. It is an interesting and important question to what extend certain 
results such as the non-exponential spin decay depend on the integrability of the model 
(3). We stress that this integrability holds very generally and does not depend on a specific 
choice of the hyperfine coupling constants Ai. 

4.2 Different types of initial states 

In the numerical simulations to be described below the electron spin is initially in a single 
tensor product state with the nuclear spin system, 

\^{t = 0)) = |Ve() ® li^nuc) , (21) 

i.e. the electron spin described by Itpei) is initially uncorrelated with the nuclear spins. 
However, there is still quite a variety of possibilities for the initial nuclear spin state 
\tpnuc)- A simple choice would be just a tensor product of eigenstates with respect to a 
given quantization axis, say, the 2;-direction, 

li^nuc) = I T)l I i)2 I i)3 ® • • • I T)iV , (22) 

where we have for simplicity assumed the nuclear spins to be of length 1/2. A nuclear 
spin state close to the above form can be generated experimentally by cooling down the 
nuclear spins in a strong external magnetic field. The strong magnetic field provides a 
quantization axis and suppresses dipolar interactions changing the spin projection along 
the field axix. Then, due to spin-lattice relaxation processes, the nuclear spin system will 
end up in a state of the type (22). 

A more general possibility of a nuclear spin state is a tensor product state but with 
individual polarization direction for each spin, 

IV'n«c) = (ail T)i + h\ i)i) • • • (ajvl T)jv + «>jv| 1)n) , (23) 

where the complex numbers a^, bi parametrize the spin state and are subject to the nor- 
malization condition |a,|^ -|- = 1. A yet more general state of the nuclear spins is a 
superposition of tensor product states, 

|^„„c) = ^aT|T), (24) 

T 

where the summation runs over all tensor product states of the form (22), i.e. over a 
complete basis of the underlying Hilbert space. If one works in a subspace corresponding 
to a given value of J^, this summation has to be accordingly restricted. For a nontrivial 
choice of the amplitudes ay states of the above form can in general not be expressed 
as tensor product states, whatever basis one would choose in the space of each nuclear 
spin. Therefore, such states are in general correlated, or, using the language of quantum 
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information theory, entangled [63, 64]. A particular class of correlated nuclear states is 
obtained when the complex amplitudes ay are chosen at random, only restricted by the 
normalization condition |q!tP = 1- In the following we will refer to this type of states 
as randomly correlated states. As wc shall see, the type of initial state of the nuclear spin 
system has a profound impact on the electron spin dynamics. 



4.3 Numerical results for electron spin dynamics: product states 
versus randomly correlated states 

Here we review our reccmt numerical studies [32] of electron spin dynamics in quantum 
dots modeled by the Hamiltonian (3). These investigations are based on exact numerical 
diagonalizations whose technical details we summarize below. 



4.3.1 Numerical method and modeling 

Our simulations of the electron spin dynamics arc performed by numerically diagonalizing 
the Hamiltonian (3). The numerical diagonalization makes use of the fact that the projec- 
tion of the total spin J = S + /j on the direction of the external field (usually chosen 
as the z-axis) is a conserved quantity leading to a block-diagonal structure of the Hamil- 
tonian matrix. Therefore it is convenient to work in subspaces of a given value of . The 
Hamiltonian is then diagonalized within such a subspace, and the time evolution of a given 
initial state is obtained from the eigensystem data. For a numerically exact simiilation of 
the time evolution of a general initial state one generally needs the full eigensystem, i.e. 
all eigenvalues and corresponding eigenvectors. If the initial state involves more than one 
of the above invariant subspaces of the Hamiltonian the time evolution in the different 
subspaces can be superimposed. We note that this method for solving for the quantum 
mechanical time evolution is non-iterative and can therefore be extended to very large 
times. On the other hand it requires the full eigensystem in a given invariant subspace 
of the Hamiltonian, and it is the dimensions of these subspaces that limit our numerical 
investigations. To reduce the numerical demands we will consider in the following nuclear 
spin of length I = \- The dimensions of the invariant subspaces increase with decreasing 

starting from the maximum value = {N -\-l)/2 according to a binomial distribution 
where N is the number of nuclei considered. The dimensions of the invariant subspaces 
become largest for the minimum value of |J^| G {0,1/2}. For the latter case we could 
simulate the full time evolution for systems with up to fourteen nuclear spins. 

In our simulations we will assume a spherical quantum dot using the following specific 
modeling. A given number N of nuclear spins is contained in a sphere of radius 

where no = 1/wo is the density of nuclei. The electron wave function is given by 



where the parameter a describes the confinement of the electron in the dot due to an essen- 
tially harmonic potential (with possibly small anharmonic corrections). In the following 
we shall use a = 2 such that the electron is reasonably confined in the sphere of radius 
R. The question of different types of confining potentials was investigated in Ref. [30]. 
Moreover, we shall use the material parameters of GaAs with no ~ 45.6nm~^. Therefore a 
typical quantum dot contains about N = 10^ nuclei. To mimic their spherical distribution 
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also in systems of smaller size used in our simulations we choose the radial coordinate 
of the i-th nucleus according to 



with i ranging from 1 to N. The results to be presented below are obtained for an 
antiferromagnetic sign of the hyperfine coupling, ^ > 0. 

4.3.2 Results for electron spin dynamics 

Fig. 4 shows results for a system of A/' = 14 nuclear spins. This is the largest system 
size for which we have been able to treat the electron spin dynamics in the presence of 
an unpolarized nuclear system. The upper left panel shows the expectation value {S^{t)) 
as a function of time for an initially fully polarized nuclear system with the electron spin 
pointing opposite to it in the negative z-direction. In the following panels the polarization 
of the nuclear system is successively reduced by lowering the value of in the initial 
state. The case of a fully unpolarized nuclear spin system is reached in the bottom right 
panel with = —1/2. Since the value of the z-component of the total spin J is fixed the 
expectation values of the transversal components (S^) and {S'^) vanish. In all simulations 
shown in Fig. 4 the electron spin is initially in a tensor product state with the nuclear sys- 
tem. The nuclear spins themselves are initially in a randomly correlated state as described 
in section 4.2. 

In all cases, = decreases in magnitude. With decreasing polarization 

the decay becomes more pronounced, and the oscillations accompanying this process get 
suppressed. Note that it is the decay of the envelope in these graphs but not the fast 
oscillations themselves that signals the decay of the spin. The distance between two 
neighboring maxima of the oscillations can depend slightly on the initial state and the 
coupling constants in the Hamiltonian. However, a good estimate for this effective period 
is usually given by T = 47r?i/^ since j4/2 is an estimate (neglecting quantum fluctuations) 
for the width of the spectrum, i.e. the difference between the largest and the smallest 
eigenvalue of the Hamiltonian. 

When the nuclear spin system is initially in a randomly correlated state the time 
evolution of {S^(t)) is very reproducible in the sense that it depends only very weakly on 
the particular representation of the initial random state. This is illustrated in Figs. 5 and 6 
where the results of different randomly correlated initial nuclear spin states are compared 
for two different system sizes and degrees of polarization. 

This behavior of randomly correlated initial states sharply contrasts with the time 
evolution of simple tensor product nuclear spin state. Fig. 7 shows the time evolution of 
the electron spin for two initial tensor product states for the same system size and degree of 
polarization as in Fig. 6. The comparison of these two figures demonstrates the significant 
difference in the electron spin dynamics for these two types of initial conditions in the 
nuclear system. In the case of tensor product initial states the time evolution depends 
significantly on the concrete initial condition, and the decay of the electron spin occurs 
typically more slowly than in the case of an initially randomly correlated nuclear spin 
system. 

4.3.3 Spin decay and quantum parallelism 

In the two left panels of Fig. 8 we show the time evolution of the electron spin averaged over 
all nuclear tensor product states for two different system sizes and initial polarizations of 
the nuclear spins. The two right panels show the corresponding data for a single randomly 
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(27) 
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correlated initial nuclear state. Comparing those plots one sees that this data is very close 
to the time evolution of a randomly correlated state. This observation is also made for other 
system sizes and degrees of polarization and constitutes an example of quantum parallelism 
[64]: The time evolution of each initially uncorrelated (and therefore c;lassical-like) nuclear 
state is present in the evolution of a linear superposition of all such states. In other 
words, the time evolutions of all uncorrelated classical-like states are performed in parallel 
in the time evolution of the randomly correlated state. An experimental consequence of 
this observation is that if the electron spin dynamics would be detected on an array of 
independent quantum dots one could not distinguish whether the nuclear spin system in 
each dot was initially randomly correlated or in an uncorrelated tensor product state. In 
other words, the spin dynamics of a randomly correlated pure state of the nuclear system 
in a single dot cannot be distinguished from a mixed state of an ensemble of dots. 

As seen above, for randomly correlated initial nuclear states the time evolution of 
the electron spin does practically not depend on the concrete realization of the random 
nuclear state and mimics closely the average over all tensor product initial conditions. 
This observation relies on the cancellation of off-diagonal terms a^ax' {il',T\S{t)\ -(J-, T"), 
T 7^ T', due to the randomness in the phases of the coefficients ax- In this sense our system 
has a self-averaging property. This can be checked explicitly by reducing this randomness. 
The left panel of Fig. 9 shows the time evolution of a randomly correlated state where the 
amplitudes ax are restricted to have a non- negative real and imaginary part. This time 
evolution turns out to be similarly reproducible as before, i.e. it does not depend on the 
concrete realization of the initial random state, but it is clearly different from the former 
case since the cancellation of off-diagonal contributions is inhibited. For comparison we 
show in the right panel of Fig. 9 data where the amplitudes in the initial nuclear spin 
state have a random phase but are restricted to have the same modulus. Here the proper 
averaging process takes place again. 

The results described so far were obtained in certain subspaces of and for the form 
of coupling constants Aj as induced by the quantum dot geometry. However, our findings 
do not depend on these choices. We have also performed simulations were the initial 
state has overlap in the full Hilbert space. For a randomly correlated initial nuclear spin 
state the only difference is that now also transverse components {S^{t)), {S'^{t)) of the 
electron spin evolve. However, these are tiny in magnitude and oscillate around zero. 
For an initial tensor product states these transverse components can become sizable, and 
the time evolution again strongly depends on the concrete initial tensor product state. 
Moreover, as mentioned earlier, the exact form of the coupling constants is also not crucial 
as long as they are sufficiently inhomogeneous. For instance, we obtain qualitatively the 
same results if we choose the coupling parameters randomly from a uniform distribution. 

4.4 Decoherence and the generation of entanglement 

In circumstances of quantum information processing the decay of a qubit is usually viewed 
as some 'decoherence' process due to the environment attacking the quantum information. 
As seen above, the spin decay is generically slower if the spin environment is initially in a 
uncorrelated state. This finding suggests that it is advantageous for protecting quantum 
information to disentangle the environment that unavoidably interacts with the qubit 
system. 

A 'decoherence' process of the above kind can be viewed as the generation of entan- 
glement between a qubit and its environment. The system investigated here provides an 
illustrative example for this statement. The entanglement in the total state \ijj{t)) between 
the central electron spin and its environment can be measured by the von-Neumann en- 
tropy of the partial density matrix where either the electron or the environment has been 
traced out from the pure-state density matrix [65]. Tracing out the nuclear 
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system we have 



n,m-f h + iS^t)) {S+{t)) 



This matrix has eigenvalues A± = 1/2 ± and the measure of entanglement reads 

E{\ijj{t))) = — A+ log A+ — A- log A_. Thus, the formation of expectation values ^ 
1/2 (or, in the case of fixed J^, just \{S^{t))\ ^ 1/2), is a manifestation of the entanglement 
between the electron spin and the nuclear spin system. The maximum entanglement, 
E = log 2, is achieved if the electron spin has decayed completely as measured by the 
expectation values of its components, {S{t)) =0. The generation of quantum entanglement 
between the electron spin and the nuclear spin system signaled by a reduced value of {S{t)) 
is a main and crucial difference between the quantum system studied here and its classical 
'counterpart' described by a system of Landau-Lifshitz equations. These equations can 
be obtained from the Hciscnbcrg equations of motion for the quantum system, dS/dt = 
i[H, S]/h, dli/dt = i[H, Ii]/h, by performing expectation values of both hand sides within 
spin-coherent states and assuming that the expectation values of all operator products 
factorizes to products of expectation values. This procedure becomes exact in the classical 
limit [66] . The resulting equations do not contain operators any more but just describe the 
dynamics of three-component vectors (classical spins) of fixed length. We have performed 
simulations of such a classical spin system by solving the Landau-Lifshitz equation via 
the fourth-order Runge-Kutta scheme. As a result, the central classical spin performs an 
irregular chaotic motion which does not show any similarity to the results for the quantum 
spin- 5 case. In particular all qualitative features of quantum effects such as the generation 
of entanglement (signaled by a decay of spins as measured by their expectation values) 
are not present in such a time evolution. Therefore the Landau-Lifshitz equation provides 
only a rather poor description of the underlying quantum system. 

We finally consider the nuclear spin correlator C{t) = {I^{t)P{0)), I = Y^ili, which 
can be measured directly by local NMR like measurements such as magnetic resonance 
force microscopy [67]. In a subspace of given and the electron spin pointing downwards 
initially this quantity reads C{t) = (J^ — {S^{t))){J^ -1-1/2). A realistic initial state will 
have its dominant weight in a series of subspaces with neighboring centered around 
some value. Then the time evolution of {S^{t)) is very similar in these subspaces, and the 
dynamics of the total nuclear spin can be mapped out by measuring the electron spin and 
vice versa. 



4.5 Averaging over nuclear configurations. Dephasing time for an 
ensemble of dots 

In section 4.1.1 we have seen (see also section 4.6.1) that the decay of (<S'^(t)) for the 
initial tensor product state occurs starting from the time t > hN/A, with fiN/A ~ 10~^s 
in GaAs dots. On the other hand, the electron spin precesses in the net nuclear field with 
the characteristic period ^ h\/N /A ~ 10~^-=-10~^s. Thus we see that the electron spin 
undergoes many precessions in a given nuclear field hj\[ = (for a given nuclear 

configuration) before decoherence sets in due to the non- uniform hyperfine couplings Aj. 
This behavior changes dramatically when we average over nuclear configurations [24] . Let 
us average Ct (t) = {T\S'' {t)—S^ (0) over all initial tensor product nuclear configurations 
|T). For that purpose we calculate CT{t) exactly by treating the nuclear field purely 
classically, i.e. as a c-number. Then we obtain, 

CT{t) = ^{l-cosih^t)), (29) 
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where = + ^ni. nuclear field, with /i^j^ = /i^^ + /i^^^. Again, the value 

of /ijv corresponds to a given nuclear tensor product state |T). We average Eq.(29) over 
a Gaussian distribution for hjss, i.e. over P(/iAr) cx exp(— 3/i|f/2w^). Defining Cci{t) = 
J dhNP{hN)CT{t), we obtain 

CMt) = l[l + {'^-l)e-'^''/% (30) 

Thus we get a very rapid (Gaussian) decay of Cci{t) for t which means that the 

dephasing time is VN /A. Prom the above equation we see that {S^{t)) saturates at 1/3 of 
its initial value of —1/2. An important approximation used here is the classical treatment 
of the nuclear field as a c- number (not an operator). The important qualitative point 
we want to illustrate here is the difference between the decoherence of an electron in an 
individual quantum dot compared to the case of an emsemble of dots. 

Investigations similar in spirit to the above considerations were performed by Merkulov, 
Efros, and Rosen [30] who have also considered the problem of the electron spin dephasing 
due to the hypcrfinc interaction for an ensemble of dots. As was already mentioned above, 
the field exerted on the electron spin by hyperfine interaction with surrounding nuclei 
is typically much larger a.t N ^ 1 than the field the electron provides to an individual 
nucleus. Following this observation Merkulov et al. [30] observed several time scales with 
different decay laws. At times shorter than hN/A they have considered an approximation 
where the nuclear spins are assumed to be static on the typical time scale of the electron 
spin dynamics. That is, in an individual qiia.ntuni dot the electron spin dynamics is 
approximated by a coherent rotation in the hyperfine field provided by the "frozen" nuclear 
spin configuration. Dephasing of the electron spin is then obtained by averaging over 
an ensemble of dots with individual nuclear spin configurations. Then they obtained 
essentially the same formula, see Eq.(30), with the same time scale- H^fN /A. The resulting 
dephasing times T| obtained in [30] are in reasonable agreement with experiments [14, 15]. 

Note that at later time scale {t 3> TiN /A) the authors of [30] obtained further slow 
electron spin decay by considering the variations of the nuclear field direction when the 
magnitude of the field is conserved, which is correct at N ^ 1. 

4.6 Further approximate studies of electron spin dynamics 

The exact solution for the electron spin dynamics under the hyperfine coupling (3) obtained 
in Ref. [24] and described briefly in section 4.1.1 is unfortunately restricted to the case 
of an initially fully polarized nuclear spin system. The numerical approach described in 
section 4.3 allows for a (quasi-) approximation-free treatment of systems with an arbitrary 
degree of initial nuclear polarization, but is, in particular at lower polarization, restricted 
to rather small systems. Therefore, in order to investigate the physically most interesting 
case of larger systems with initially moderate or even very low nuclear polarization, one 
needs to resort to approximations. In this section we review several recent approaches. 

4.6.1 Perturbation theory 

The contribution to the hyperfine Hamiltonian (3) coupling to the ^;-component of the 
electron spin, 

no=g^lBS'B + S'-Y,A^^ (31) 

i 

is diagonal in a basis of tensor product states with respect to the ^-axis as formulated in 
Eq. (22). The remaining part of the Hamiltonian, 

^ = E Y (S^^i + ^"^^) ' (32) 
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was treated in Refs. [24, 25] as a perturbation to Hq. Again we shall concentrate on the 
case of nuclear spins of length 1/2. Assuming the system to be initially {t = 0) in a 
tensor product state with respect to the 2;-axis (cf. (22)) with the electron spin as before 
pointing downwards, the lowest nonvanishing contribution to the electron spin dynamics 
in time-dependent perturbation theory is of second order in V. Specifically, one finds 



where the summation goes over all intermediate states \k) which are tensor product states 
of the form (22). Vik is the matrix element of V between an intermediate state \k =if 
,{...., = -1/2,...}) and the initial state \i =-IJ., {..., = +1/2,...}), and oJik = (cj - 
£k)/h, where e^, Sk are the eigenvalues of \i) and |fc), respectively, with respect to Ti.Q. 

Evaluating this lowest-order contribution for a large {N ^ 1) unpolarized system, one 
finds an universal power law for the decay of {S^{t)): For times large compared to UN/ A 
|(5^(t))| decays as fr^/"^. This is a central finding within this pcrturbativc approach and 
agrees with the perturbative limit (i.e. large Zeeman field) of the exact solution in the fully 
polarized case discussed in section 4.1.1. Note that for a weak Zeeman field Cz < A/^/N 
the part of the electron spin state which decays is of the order of the initial value, in 
contrast to the fully polarized case where this part is of order 1/A^, see Fig. 2. 

The perturbative approach has the following shortcomings: Clearly, Hq and V are 
of the same order of magnitude if one does not apply a very large; extcirnal magnetic 
field. Therefore, there is in general no small parameter justifying a perturbative expansion 
concentrating on low orders, and one would need to sum over higher orders, provided that 
this perturbative series has sufficient convergence properties. In higher orders, however, 
one encounters increasing divergences due to vanishing denominators in the perturbative 
contributions [24]. Therefore, the findings obtained from the low-order result (33) might 
appear not very reliable. Indeed, from the exact solution in the fully polarized case [24, 25] 
(cf. section 4.1.1) one finds a different (not a power) law for spin decay in the limit of low 
external magnetic field. However, there is a reasonable hope that the basic conclusion from 
the low-order perturbative approach is still correct, namely that a non-uniform hyperfine 
coupling leads to a non-exponential spin decay. 

Moreover, the perturbative approach can, for technical reasons, only deal with initial 
states where the nuclear spin system is in a tensor product state. As seen in section 4.3 
the behavior of such type of initial states depends significantly on the particular initial 
condition, in contrast to the behavior of randomly correlated initial states. 

4.6.2 Studies using Markovian approximations to the nuclear spin dy- 
namics 

Saykin, Mozyrsky, and Privman [31], and de Sousa and Das Sarma [28, 29] have performed 
investigations using, among other simplifying assumptions, Markovian approximations to 
describe the dynamics of the nuclear spins. 

In Ref. [31] the situation of an electron bound to a ^^P donor in a silicon matrix was 
studied. The electron spin interacts with the central ^^P and the surrounding ^^Si nuclear 
spins. Then the authors formulate a master equation for the reduced electron density 
matrix where the nuclear spin dynamics is governed by a Markov process. Moreover, the 
total density matrix is assumed to be separable at all times with respect to the subsystems 
given by the electron spin and the nuclear spin bath, and the density operator of the latter 
system is assumed to be time- independent. Due to the first assumption, decoherence due to 
the formation of entanglement as described in section 4.4 is excluded since the total density 
matrix is taken to be always separable [64]. As a result, the authors find an exponential 
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decay of tlic' dements of the reduced density matrix to their equihbrium values, where the 
longitudinal relaxation time Ti and the transversal dephasing time T2 (with respect to 
the direction of a weak external field) fulfill the relation Ti = T2/2. However, as already 
explained in section 4.1.1, the application of the Markov approximation in this situation 
is not justified, which might explain that this result differs qualitatively from the findings 
from the exact solution in the case of an initially fully polarized nuclear spin system (cf. 
section 4.1.1) and the perturbative result for the unpolarized case (cf. section 4.6.1), where 
a power law decay was found. 

In Refs. [28, 29] the Hamiltonian (3) was studied under the assumption of strong 
external magnetic field coupled to the electron spin and suppressing the spin-flip terms. 
Therefore, the hyperfine interaction is approximated by (31). This approach is applicable 
at an external Zeeman magnetic field which is much larger than the internal nuclear field, 
i.e. Sz ^ hujN ~ A/\/N. We recall that the part of the initial electron spin state which 
decays due to inhomogeneous hyperfine coupling is small under the above conditions, this 
part is of the order of {Tiujn/£zY < 1) see Refs. [24, 25]. The authors of Refs. [28, 29] then 
introduced the nontrivial dynamics through the dipolar interaction among the nuclear 
spins. This interaction is taken into account in a truncated version neglecting terms that 
change the total spin component in the direction of the external field, consistent with 
the earlier approximation. Thus, in this approximation the dipolar dynamics among the 
nuclear spins induce the decay of the electron spin. As a further assumption, the nuclear 
spin system is then approximated as a fluctuating fleld generated by a Markov process. 
After performing a detailed mathematical analysis of this resulting effective model [29] the 
authors conclude that the electron spin decay will occur on essentially the same time scale 
as the dipolar interaction among the nuclear spins, i.e. T„2 — 10~^ -j- 10~^s. Given the 
various assumptions leading to this effective model, this is a very natural result. 

Applying their findings to the situation of an electron bound to a ^^P donor embedded 
in an Si crystal, the authors find an interesting dependence of the spin memory time Tm (as 
measured by spin-echo experiments) on the direction of the external field relatively to the 
crystal. This effect is induced by the directional dependence of the dipolar interaction and 
provides a possibility to experimentally probe these results and the underlying assumptions 
[44]. 

4.7 Further studies and developments 

We now summarize further, mostly theoretical, work relevant to the issue of hyperfine 
interaction between an electron spin bound to a quantum dot and surrounding nuclear 
spins. 

The relaxation rate of the longitudinal electron spin component, i.e. 1/Ti, which is due 
to the interplay between hyperfine interaction and dissipative phonon processes was studied 
by Erlingsson, Nazarov, and Falko [23], and by Erlingsson and Nazarov [27]. The approach 
in Ref. [27] is a semiclassical one introducing an internal field due to the nuclear spins which 
acts on the electron spin in addition to an external magnetic field. The observed relaxation 
rate is very small, though this mechanism can prevail over the mechanism which is the 
interplay of the spin-orbit interaction and spin- independent interaction with phonons [20] , 
it happens at very low external magnetic field when the corresponding Ti time is order of 
100s. The analogous problem of the calculation of the T2 time which is due to the combined 
effect of hyperfine interaction and phonon processes was considered recently by Semenov 
and Kim [40]. These authors noticed that due to thermal fluctuations the electron can make 
spin-conserving transitions between different orbital states. Then because of either the 
energy dependence of the g-factor or different nuclear fields seen by electron spin in different 
orbital states these fiuctuations lead to electron spin decoherence (since the precession 
frequencies are different in different states). These mechanisms are only important at 
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relatively high temperatures since they are exponentially suppressed at temperatures much 
smaller than the energy distance between the neighbouring orbital states. For example, in 
the case of combined effect of hyperfine interaction and phonon processes for typical GaAs 
quantum dots the decoherence time is of order of 10s at temperature IK. 

Lyanda-Geller, Aleiner, and Altshuler [26] have investigated the dependence of the 
relaxation rate of nuclear spins in a quantum dot should on the electronic state of the 
dot and concluded that effects of Coulomb blockade should also be observable in nuclear 
relaxation in such systems. 

A proposal for using the nuclear spin system as a long-lived memory for information 
originally contained in the electron spin qubit was put forward by Taylor, Marcus, and 
Lukin [37]. In a related study, Imamoglu, Knill, Tian, and Zoller [38] have proposed an 
all-optical scheme for polarizing the nuclear spins by manipulating the electron spin. 

Saikin and Fediclikin [35] have investigated the influence of hyperfine interactions on 
gate operations within the Si:P quantum computing proposal due to Kane [5]. 

The possibility of nuclear spins forming an effective quantum dot confining electrons 
through hyperfine interaction was proposed recently by Pershin [39]. 

Theoretical studies on the decoherence of a two level systems coupled to surrounding 
spins were also presented recently by Frasca [68]. These studies stem from a somewhat 
different context but are similar in spirit the ones reported on in this article. We also 
mention a recent numerical study by Dobrovitski et al. [69] on spin dynamics stressing 
the role of entropy. There a central spin is coupled inhomogeneously to an essentially non- 
interacting spin environment where an Ising-like coupling was used. To allow for nontrivial 
dynamics the authors introduced a magnetic field perpendicular to the ^-direction of the 
Ising coupling. In a very recent paper [70] the same authors have numerically studied 
the damping of quantum oscillations in the system of two central spins. These central 
spins are coupled by strong Heisenberg exchange, this system in turn is coupled to a spin 
bath (which can be a nuclear system) through inhomogeneous Heisenberg interaction, 
a scenario smilar to the one studied in [3]. The state of the bath is initially a random 
superposition of all possible basis states. Assuming the coupling constants are random, 
in the system of A'' = 13 bath spins the authors observed a two-step decoherence process. 
Initially decoherence is very fast and after the first step the oscillations of z-components 
of the central spins Si{t) and 5*1 (t) persist and decay very slowly. It is interesting that 
the authors have managed to fit the first step fast decay by Eq.(30), see Sec. 4. 5, which 
probably confirms the self-averaging occurring for randomly correlated initial nuclear states 
described in Sec. 4. 3. 3. The second slow step of the decoherence process is presumably the 
analogue of the decay described in sections 4.1.1 and 4.6.1. 

5 Conclusions and outlook 

We have reviewed the recent literature on electron spin dynamics in semiconductor nanos- 
tructures due to hyperfine interaction with surrounding nuclear spins. This issue is of 
particular interest with respect to several proposals for quantum information processing 
in solid state systems. Although the basic Hamiltonian (3) looks rather simple it describes 
an intricate many-body problem which in general does not seem to allow for an analytical 
solution. 

In the case of an initially fully polarized nuclear spin system an exact analytical solution 
for the spin dynamics can be found. For not completely polarized nuclei approximation- 
free results can only be obtained numerically in small model systems. We have compared 
these exact results with findings from several approximation strategies such as perturba- 
tion theory and (semi-)classical approximations to the nuclear spin dynamics including 
Markovian approximations. 
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Many of those approximation arc not particularly well controlled, and they suppress 
important features of the full quantum system. For instance, the pronounced dependence 
of the dynamics of the electron spin on the type of initial condition for the nuclear sys- 
tem is not reproduced by any of the approximation strategics which have appeared so 
far. Therefore, the most obvious direction for future work is the development of possibly 
systematically controlled approximation techniques which reproduce important features of 
the full quantum dynamics and allow reliable predictions for realistic systems. A possible 
but presumably techically quite complicated route for such future work is the Bethe ansatz 
solution outlined in section 4.1.2. Such progress toward more reliable predictions for exper- 
imentally relevant situations is especially desirable because of the importance of this issue 
to several proposals for quantum information processing in a semiconductor environment, 
see in particular Refs. [3, 5], and possibly also for other scenarios in the emerging field of 
spin electronics [1, 2]. 
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Figure 1: Energy spectrum (in units of A > 0) of the Hamiltonian (3) for zero magnetic field 
as a function of the total spin quantum number J in a system of = 13 nuclear spins of 
length 1/2. For details about the modeling of coupling constants see section 4.3.1. 




Figure 2: Schematic dependence of < Sz{t) > on time t for the unpolarized tensor product 
and fully polarized nuclear states. The time scale for the onset of the decay ~ N/A is the 
same for both cases. In the fully polarized case the magnitude of the effect is 1/A^. The period 
of oscillations is of the order of ^/N j A for the unpolarized and ^ Xj A for the polarized case. 
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Table 1: Natural abundances, total nuclear spin quantum numbers /, and nuclear magnetic 
moments \ii of spin-carrying stable nuclei relevant to semiconductor materials. The values 
for ^^J are given in units of the nuclear magneton /xat = eh/2mp where rup is the proton 
mass. The data is adopted from the American Institute of Physics Handbook, Third Edition, 
McGraw-Hih, 1972. 
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Figure 3: The dependence of |ap averaged over time (< >) on the external Zeeman 

field €z for a fully polarized nuclear state. The resonance occurs at le^j = A/2, and the width 
of the resonance is ~ A/^/N, which is much smaller than the initial gap A/2. 
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Figure 4: The time evolution of the electron spin in a system of = 14 nuclear spins of 
length 1/2 for different degrees of polarization of the randomly correlated nuclear system. The 
hyperfine coupling constants are induced by the quantum dot geometry. In all simulations the 
electron spin is initially pointing downward in a tensor product with the nuclear system. In 
the top left panel the nuclear spins are fully polarized in the initial state with the electron spin 
pointing opposite to them (J^ = 13/2). In the following panels the number of flipped nuclear 
spins in the initial state is gradually increased. The case of an initially fully unpolarized (but 
randomly correlated) nuclear system is reached in the bottom right panel (J^ = —1/2). Here 
and in the following we take spins to be dimensionless, i.e. measured in units of h. 





Figure 5: Electron spin dynamics for an initially randomly correlated nuclear spin system 
for N = 19 nuclear spins with a moderate degree of polarization (J^ = 6). The data of six 
completely independent random realizations of the initial nuclear spin state are shown. The 
resulting electron spin dynamics is practically independent of the realization of the initial state 
for this type of initial condition. 
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Figure 6: Data of the same type as in Fig. 5 for = 29 nuclei and = 12. Again, the electron 
spin dynamics is practically independent of the realization of the random initial nuclear state. 
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Figure 7: Electron spin dynamics for = 29 nuclei and = 12. Here the initial state of 
the nuclear spins is given by individual tensor product states. Different initial tensor product 
states clearly lead to a significantly different time evolution of the electron spin. As seen in 
Fig. 6 this is strikingly different from randomly correlated initial conditions. 
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Figure 8: The left two panels show the electron spin dynamics averaged over all possible initial 
tensor product states for two different system sizes and degrees of polarization. The right two 
panels show the corresponding data for a single randomly correlated initial condition for the 
nuclear system. As the comparison shows, the time evolution for the randomly correlated 
nuclear spin system closely mimics the average over all tensor product initial states. 
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Figure 9: Time evolution of {S^{t)) for two types of initially randomly correlated nuclear 
spin states. In the left panel the amplitudes are restricted to have non-negative real and 
imaginary part, while in the right panel they have all the same modulus but completely random 
phases. 
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